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We derive the exact partition function for a discrete model of random trees embedded in 
a one- dimensional space. These trees have vertices labeled by integers representing their 
position in the target space, with the SOS constraint that adjacent vertices have labels 
differing by ±1. A non-trivial partition function is obtained whenever the target space 
is bounded by walls. We concentrate on the two cases where the target space is (i) the 
half-line bounded by a wall at the origin or (ii) a segment bounded by two walls at a finite 
distance. The general solution has a soliton-like structure involving elliptic functions. We 
derive the corresponding continuum scaling limit which takes the remarkable form of the 
Weierstrass p function with constrained periods. These results are used to analyze the 
probability for an evolving population spreading in one dimension to attain the boundary 
of a given domain with the geometry of the targets (i) or (ii). They also translate, via 
suitable bijections, into generating functions for bounded planar graphs. 
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1. Introduction 




Fig.l: A sample rooted labeled planar tree, with root (top) vertex labeled 
by 2. Neighboring vertices have labels differing by ±1. These labels may be 
viewed as positions on a target integer line as indicated. 

In this paper, we consider a simple model describing the embedding in one dimension of 
a random tree. More precisely, we consider random rooted trees whose vertices are labeled 
by integers representing their possible discrete positions in a one-dimensional target space. 
We moreover impose that two neighboring vertices on the tree have labels differing by 
+1 or —1, which allows to view the edges of the tree as rigid segments of unit length 
embedded in the real line (see Fig.l for an illustration). We choose to consider the case 
of so-called planar trees, i.e. we count as distinct all trees obtained by permuting any 
two descendent subtrees at a given vertex. This model is nothing but a discrete version 
of the so called one-dimensional Brownian snake which is used to describe branching 
processes, for instance the spreading of a population in a one-dimensional target space. In 
this language, we may think of our rooted trees as representing "genealogical trees" for the 
lineage of an initial individual (materialized by the root vertex), while the labels represent 
the positions in space of all the descendent individuals. The spreading process is modeled 
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here by demanding that each individual lives at distance one from its parentS. We will 
discuss this interpretation in detail in section 5 below. 

Alternatively, we may view this model as a statistical Solid-On-Solid (SOS) model with 
heights given by the labels, and whose base space is a random tree. The SOS rule of having 
neighboring heights differing by ±1 is responsible for the "integrability" of the model. 
However, as opposed to the two-dimensional lattice case where a roughening transition 
takes place, we expect here that the discrete nature of the heights is eventually irrelevant 
when working with large tree-like base spaces. 

We shall consider a statistical sum over all such trees with a weight g per edge, and 
with a fixed position, say n, of the root vertex. With no bounds on the labels, the partition 
function is trivial as it amounts to counting rooted planar trees with N edges (in number 
ctv where cn = ( 2 ^) / (N + 1) is the iV-th Catalan number), each of which gives rise to 2^ 
possible embeddings. Similarly, the width wn = a/ (n 2 ) n for the fluctuations of the labels 
n in a random tree of size iV and with root label is easily obtained as 

2 _ i ( 4" . 

hence wn ~ (nN/A) 1 / 4 for large N. 

The problem becomes more interesting in the presence of a wall, say at position — 1, 
which amounts to imposing that all labels be non-negative. Similar trees were introduced 
under the name of "well-labeled trees" in Ref. [|J in connection with the enumeration of 
rooted planar tetravalent graphs, with vertex labels representing the (necessarily non- 
negative) geodesic distance on the graph to the root vertex!!. The explicit form of the 
corresponding partition function as a function of g and n was given in Ref.fH] in the 
context of planar graph enumeration and will be recalled in the next section. A remarkable 
outcome of this solution is the emergence of discrete soliton-like expressions, suggesting an 
underlying integrable structure. 

The purpose of this paper is to study the statistics of trees with labels now belonging 
to a finite set, say {0, 1, • • • , L}, which amounts to having two walls at positions —1 and 

4 The particular discrete spreading rule chosen here should not affect the universality of the 
continuum answer. 

5 In this reference, a slightly different constraint is imposed on the labels, demanding that 
neighboring vertices have labels differing by 0,±1. Our results will be easily extended to this 
modified case in Sect. 6 below, with no fundamental difference. 
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L + 1 in the embedding target space. In the language of evolving populations, introducing 
walls gives access to the probability for the population to attain pre-defined boundaries 
or to remain confined within a pre-defined connected domain. The two-wall situation 
corresponds to the generic case where this domain is compact, i.e. is a segment (with two 
boundaries). Such boundary conditions correspond to the so-called Restricted Solid- On- 
Solid (RSOS) version of the problem, in which heights are restricted to belong to a finite 
segment. 

The paper is organized as follows. In Sect. 2, we introduce a master equation for the 
partition function of labeled rooted trees and recall its solutions both without walls and 
with one wall. Sect. 3 is devoted to the derivation of the two-wall solution, involving elliptic 
functions. The corresponding continuum limit is derived in Sect. 4 and expressed in terms 
of the Weierstrass p function with constrained periods. As an application of these results 
we study in Sect. 5 a particular stochastic process describing the evolution of a population 
which spreads in one dimension. We discuss in Sect. 6 the solution of a slightly different 
problem corresponding to a dilute SOS version in which neighboring vertices of the tree 
have labels differing by ±1 or 0. We gather a few concluding remarks in Sect. 7, where 
we discuss the integrability of our models in particular in connection with planar graph 
enumeration and matrix models. The precise connection between labeled trees and planar 
graphs is further detailed in Appendix A. 



2. Enumeration of labeled rooted trees 

Let R n denote the generating function for labeled rooted planar trees with a root at 
position n. By a decomposition according to the possible local environments of the root 
vertex, characterized by the sequence of labels n — 1 or n + 1 of its adjacent vertices (see 
Fig. 2), we immediately get the equation 

1 — d(Rn+l + Rn-l) ^ 

Note also that from their combinatorial definition as counting functions, the i? n 's are 
required to have a series expansion in powers of g starting as R n = 1 + O(g), a condition 
sufficient to determine all the -R n 's from eq. (|2.1| ) . The relation (|2.1| ) is valid for all accessible 
values of n and it may be supplemented by boundary conditions to account for the possible 
presence of walls. We choose those conditions among three categories: no wall, one wall and 
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Fig. 2: A sketch of the master equation (2.1). A tree contributing to R n 
is decomposed according to the sequence of descendents of its root (labeled 
by n) . Each descendent vertex is arbitrarily labeled by n ± 1 henceforth the 
descendent subtrees are generated by R n ±\ accordingly. Each edge connected 
to the root is given a weight g. The summation over all possible configurations 
produces the r.h.s. of eq.(2.1). 

two walls, corresponding to restrictions on the allowed labels n, respectively no restriction, 
n > and < n < L. 

In the absence of wall, all the R^s are equal due to translational invariance, to a 
function R satisfying R = 1/(1 — 2gR) and R = 1 + O(g), namely: 

R = R(g) - '"f^ = £ 9 N 2 N c N (2.2) 

^ N=0 

This formula displays the critical value g c = 1/8 of g, while the coefficient of g N clearly 
counts rooted planar trees with N edges (cat) with 2 N possible embeddings. 

In the presence of one wall at position —1, we must write R-i = and consider the 
equation (|2.1| ) only for n > 0. This system may be solved order by order in g using as 
a seed the vanishing of all the coefficients of the series for and the order zero values 
R n = 1 + O(g) for all n > 0. In a more global way, the solution was worked out in 
Ref.0 by replacing the condition of existence of a power series expansion for each R n 
by the condition that R n — > R at large n with R as above. Indeed, it is clear from the 
combinatorial definition that R — R n = 0(g n+1 ) as the wall at position —1 may not be 
reached with less than (n + 1) edges from position n. The solution of eq.( |2.1| ) with these 
boundary conditions reads: 

R n = R Un n+4 , u n =x 11 ^ 1 -x~ 11 ^ 1 (2.3) 

Un+lUn+3 



4 



with R = R(g) as in eq^^) and where x is the solution of 

x + - + 2 = -L (2.4) 
x gR 

with, say, modulus less than 1, namely x = x(g) = (1 — (1 — 8g)~z)/(l + (1 — 8g)%). Note 
that x is real for all g < 1/8 and admits a convergent series expansion in g with positive 
integer coefficients. The small g behavior x(g) = g + 0(g 2 ) ensures the above property that 
R n = R + 0(g n+1 ). The solution ( |2~3l ) is readily checked by noting that eq.(|0]) translates 
into the following trilinear equation for the it's 

U n U n+2 U n+ 4 = —U n+1 U n+2 U n+3 + ^i?(w n _iM n+ 3« n+ 4 + U n U n+1 U n+5 ) (2.5) 

it 

easily verified upon substituting gR = x/(l + x) 2 and 1/R = (1 + x 2 )/(l + x) 2 . The 
particular form of the solution ( |2~3| ) was identified in Ref. j3| as a stationary one-soliton 
solution to the KP equation Q. 

In the presence of two walls, we must write R-i = Rl+i = and consider eq. ( |P| ) 
for < n < L. This system may again be solved order by order in g from the boundary 
conditions that all coefficients of R-± and Rl+i vanish and R n = 1 + 0(g) for all n, 
< n < L. For a finite value of L, one may also eliminate all but one R n from the finite 
set of algebraic equations ( |2.1|) so as to obtain an algebraic equation of degree [(L + 3)/2] 
for each R n , with a unique solution such that R n = 1 + O(g). In a more global way, 
we intend to generalize the solution ( |2.3|) to this two-wall case. This is done in the next 
section. 



3. Two-wall solution 

3.1. General solution via elliptic functions 

From the one-soliton structure of the solution (12. 31), it is natural to look for a more 



general soliton-like solution to the equation ( |2.1|) in the "elliptic" form 



Rn — R- 

Un = (x^ - x~^) JJ(1 - q j x n+1 ){l - 



oo j (3.1) 



3=1 



with an additional free real parameter q such that \q\ < 1. For q = 0, we recover the 
solution (|2.3|) provided R and x are given by eqns. (|2.2| ) and (|2.4j), depending on g only. As 



we shall now see, for a general q, we may tune R and x in eq.( |3.1|) as functions of both g 
and q so as to satisfy the equation ( J2.1| ). The above solution clearly satisfies R-i—0 and 
the value of the free parameter q may finally be adjusted so as to ensure Rl+i = 0. 

To fix the functions R and x in terms of g and q, we again write eq. (|2.1|) in the form 
of eq . ( p75|) , and note that from the definition of u n , we have it_i = and U-2-k = —Uk- 
Hence, taking n = — 1 in eq.(|2.5|), we get gR 2 = u\/u% while taking n = —2, we obtain 
gR = (uq/ui) 2 , which generalizes eq. (|2.4| ). This leads to: 

9 - ^ (3-2) 

which implicitly determines x(g, q) as a function of g and q. More precisely, for a fixed 
q with < q < 1, this equation has four solutions for 0<(7<l/8,a real positive one 
x o(9i q) with, say, modulus less than 1 together with its inverse X2(g, q) = l/xo(g, q), and 
a solution on the unit circle x±(g, q) with, say, positive imaginary part together with its 
inverse xs(g, q) = l/xi(g, q). For a fixed q with — 1 < q < 0, we have the same pattern of 
solutions provided g > 1/8 and does not exceed some upper bound. When g = 1/8, all 
solutions coalesce to x = 1 independently of q. For a definite value of the position L of 
the second wall, the proper choice of solution will be discussed below. In particular, the 
presence of two walls at a finite distance increases the radius of convergence of the -R n 's 
as series of g as it reduces the entropy of configurations. This requires exploring values of 
g > 1/8. Note that for q = 0, we have xo(g, 0) = x(g), corresponding to the solution of 
section 2, while x x (g, 0) = (1 + i(l - 8#)*)/(l - i(l - 8g)%). 
The function R is now given by 

R=^P- (3.3) 

with two determinations R(°\g,q) and R^(g,q) corresponding respectively to the sub- 
stitutions x = xo(g,q) and x = x\(g,q) in the w n 's, or equivalently X2(g,q) and xs(g,q) 
respectively as the change x — > 1/x amounts to u n — > —u n , leaving R and all i? n 's invari- 
ant. It remains to verify that, for these particular choices of x and R, the w n 's of eq.( |3.1| ) 
actually satisfy the identity (|2.5|) for all n. Let us introduce the notations q = e 2l7rr , 
x n+1 = e 2l ' KZ so that u n = 0\(z) where Q\ is the (unnormalized) Jacobi theta function with 
nome q and argument z: 

6 x (z) =0(z\t) = 2zsin(7rz) JJ(1 - 2q J cos(2tvz) + q 2j ) (3.4) 

j>i 
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We also introduce the notation x = e 2l7ra , so that eq. (|2.5|) , together with the particular 
choices ( |3.2|) and (|3.3|), translates into a theta function identity 

2 

&1 ^ ' 1 {z)6 1 (z + 2a)6 1 {z + 4a) = 6 ^\o 1 {z + a)6i{z + 2a)d 1 {z + 3a) 



Oi(a) J v ' iV ' v ' 0i (2a) 

+ #i(z - a)6>i(z + 3a)6>i(z + 4a) + 0i(z)0i(z + a)Q\{z + 5a) 

(3.5) 

To prove it, we note that both hands have the same transformations when z — > z + 1 and 
z — > z + r, due to the properties 9\{z + 1) = — 6>i(z) and 6>i(z + r) = — g~5e~ 2l7rz 1 (z). 
Taking the ratio of rhs/lhs of eq. (|3.5| ), we get an elliptic function, with poles possibly at 
z = 0, —2a, —4a in a fundamental cell. One easily checks by examining the rhs at these 
values that the corresponding residues all vanish and therefore the ratio is a constant, 
easily shown to be 1 by taking for instance its value at z = —a, which completes the proof 
of the identity. With these notations, we have the relations 

9 = 0S(2«) (3 - 6) 

and 
while 

= 0i(2a) 3 1 ((n + l)a)0 1 ((n + 5)q) 
n 0i(a) 2 0!(4a) 9 1 {{n + 2)a)e 1 {{n + A)a) { ' 



with 6*i (z) as in eq.(pT4]). As this solution involves elliptic functions, it is natural to study 
its transformation under the modular transformation r — > — 1/r. From the transformation 
9\{z/t\ —1/r) = — z(— zr) 1 / 2 e l7r2 / T 0i(z|r), we find that the physical quantities such as # as 
given by eq.( |3.6| ) above and R n as given by eq. (p78|) are invariant under (a, r) — > (a, f) = 
(a/r, —1/r) or equivalently (x, g) — > (x, g) with x = x 1 ^ and g = e _2l7r / T . This is not the 
case for intermediate factors such as R or any of the «„'s taken independently. We deduce 
from these properties that the same solution is reached by the parameters (x, q) and by 
their modular transforms (x,g). 

It is interesting to study the modular transformation of the solutions Xi(g,q), i = 
0, 1, 2, 3 of eq.(|3.(3|) as introduced above. Clearly, for a fixed g, the modular invariance of 
the r.h.s. of eq. (|3.6| ) implies that xt(g,q) = x^(g, g) 1//r is a valid solution when g — > g, 
hence we may write x^(g, g) 1//T = x ff u)(g,q) for some permutation o G S4. Iterating this 
transformation, and using q = q while rf = —1, we deduce that x a 2^(g,q) = Xi(g, g) _1 , 



which brings us back to the same solution we started from, but with the determination 
Xi(g, instead of Xi(g, q). Therefore a is a circular permutation of the four indices. 

For q < 0, the modular transform q becomes complex, and we will make no use of 
the above remarks. On the other hand, for < q < 1, the modular parameter r = it may 
be taken purely imaginary, and the modular transformation reduces to t — > 1/t hence q is 
also real and in the range (0, 1). It is then easy to see that the modular transformation 
sends the solution xo(g, q) to x\(g, q) and xi(g, q) to X2(g, q) which leads to the same -R n 's 
as xo(g, q). In other words, the same physical solution may be described by either some q 
and the real determination xo(g, q) or by its modular transform q and the determination 
on the unit circle xi(g,q). 

3.2. Boundary condition at n = L + 1 

We now implement the two-wall boundary condition described above, in which we 
require Rl+i = 0, or equivalently ul+5 = 0. This is achieved by demanding that x L+6 = 
q m for some m G Z. From the above study, we have at our disposal absolutions either 
real positive and smaller than 1, or on the unit circle with positive imaginary part (we 
discard the equivalent l/x-solution). This leaves us for positive q with two possibilities: 
(i) x = e 2l7rx +e, k = 1,2,... corresponding to m = and (ii) x = m = 1,2,... 

while at negative q, only the solution (i) survives. As the solution R n must be positive 
by definition for n = 0, 1, L, it is easily verified that only k = 1 in case (i) and m = 1 
in case (ii) are admissible to prevent sign changes for the « n 's and the -R n 's. In a more 
physical language, higher values of k or m correspond to higher modes with oscillations in 
the range [0, L]. Taking for instance m = 2 corresponds to a first vanishing of R n at the 
coordinate n = L/2 — 2. 

For any given L, we may always pick the solution (i) and define qL(g) as the unique 
real solution of 

x 1 (g,q L (g)) = e^ (3.9) 
Introducing the notation L(g) for the value of L such that xi(g, q = 0) = e L + 6 , namely 

L(g) = -A -r- - 6 (3.10) 

arctan ( (1 — 8g) 4 j 

we have that qh{g) > for L > L(g) and qL{g) < for L < L(g), also valid for g > 1/8 
with the convention that L(g) = 00 in this range. The length L(g) may be taken as a 
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measure of the typical extent in the embedding space of the random trees at a fixed value 
of g, in the absence of walls. 

On physical grounds, we expect a qualitative change of behavior to occur at wall 
distances L of the order of L(g). For L » L(g), the tree behaves as a compact object of 
typical extent L(g) "diffusing" between the walls and feeling them only when approaching 
at distances of the order of L(g) and smaller. For L << L(g), the walls strongly squeeze 
the tree and L is the only relevant scale of the problem. In the first regime, the profile 
{R n }o<n<L, always maximal in the middle, is mainly constant and decreased significantly 
only at distances of order L(g) from the walls. In the second regime, the profile will vary 
over the whole range between the walls. Alternatively, for a fixed value of L, we may invert 
these conditions into qi,(g) > for g < gi, and qh{g) < for g > gi,, where 

9L = Ul-t^ (-?-)) (3.11) 



.£ + 6. 

Note that < g L < 1/8 for all L = 0,1, 2, ... 

Beside the solution ( |3.9| ) above, we have another possibility for the choice of q namely 
q = q L (g) which solves the condition (ii) 

xo(g,q' L (g)) = q' L (g)^ (3.12) 
This alternative solution exists for all L > L{g) and is a priori distinct from the previous 



solution (|3.9| ). In practice, the physical solution for fixed g and L is unique as one easily 
checks that the solution ( |3.9| ) is the modular transform of ( |3.12j ), namely qL,(g) = q'lG?) 
and x x {g, q L (g)) = x Q (g, q'^g)) 7 ^ with q' L (g) = e 2l7TT '^ 9 l The two choices Q or ( gUp 
therefore lead to the same physical quantities R n . 

To obtain the combinatorial series expansions for the -R n 's in g it is simpler to work 
with the solution (|3.12| ). This is always possible as g is small (hence we may work in the 
regime g < gL). Moreover, we have xo(g, q' L (g)) = g + 0(g 2 ) while q' L (g) = g L+G (l + 0(g)). 
This shows in particular that the present R n and the former R(g) of eq.( |2.2| ) have the same 
expansion up to order Min(n + 1, L — n + 1) in g, as expected. 

On the other hand, to have a more global approach it is best to work with the solution 
which is valid for any g. We may then use the relations ( |3.6| )-( |3T8| ) with a = l/(L + 6) 
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to parametrize the solution with q as follows 
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with #i as in eq. fl3.4| ). This form displays clearly the symmetry R n = Rl-u expected from 
the symmetry of the problem, as a consequence of the relation #i(l — z) = 6\{z). 
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Fig. 3: Plot of g(q) as given by eq.( [3-13|) for L = 6. The function decreases 
from a maximum value g c (L) at some negative value q c (L) down to at q = 1. 
At q = 0, we have g = gL as in eq.( |3.11| ). We have gL < | < g c (L). 



For fixed L, the function starts from g = at q = 1 and increases as g decreases, 
passing through of eq.( |3.11| ) at g = 0, and reaching a maximum value g c (L) > 1/8 
attained at some negative value of q = q c (L) (see Fig. 3 for illustration). The quantity 
g c {L) is nothing but the radius of convergence of all the series in g appearing in the 
problem, and governs the leading growth of the number of configurations as a function 
of the number iV of edges in the tree as g c (L)~ N . We have for instance g c (^) = 1/4, 
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Fig. 4: Typical profiles {R n }o<n<L for L = 100 as obtained from eq.( |3.13| ) 
for three particular values of q: (a) a positive value realizing g < where the 
profile is flat except for a region of extent L(g) from the walls (here q = .25 
and L(g) ~ 18) (b) <? = (# = ^l) and (c) q = q c (L) (g = g c (L)) where the 
profile is slightly peaked in the middle. 

g c (2) = 3 - 2y/2, g c (3) = 4/27. We also see that g c (L) -^1/8 when L -> oo. The branch 
of the solution for g < g c (L) is discarded as unphysical. 

To conclude this section, let us use the solution ( |3.13| ) to display for fixed L the exact 
profile {R n }o<n<L for some particular values of q, namely 

(a) a positive value of q realizing g < g^, in which case the profile is flat except for a 
region distant by typically L(g) from the walls 

(b) the value q = where g = gi and all theta functions degenerate into trigonometric 
functions 

(c) the negative value q = q c (L) corresponding to g = g c (L), where the profile varies 
over the whole range [0, L] and is slightly peaked in the middle 

These profiles are represented in Fig. 4 for L — 100. 
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4. Continuum limit 

For starters, let us derive the continuum limit of the one- wall solution ( |2.3| ). It is 
reached by letting g — > 1/8 and n — > oo simultaneously as: 

9=l(l-e 4 ), n=- e (4.1) 

with e — > playing the role of the inverse of a correlation length. This leads to 

x(9) = = ^ (4.2) 

and finally, expanding R n up to order 2 in e: 

# n = 2(l-e 2 W(r)) (4.3) 

where U is the scaling function 

W(r) = 1 + — |— (4.4) 
Sinn (rj 

describing the "repulsive potential" felt by the tree as a function of the rescaled distance 
r from the wall at r = 0. 

In the presence of the second wall, we let in addition L — > oo by keeping the quantity 
uj = (L + 6)e/2 fixed. This ratio of the two characteristic lengths of the problem, namely 
L and the typical extent of the unconstrained tree L(g) ~ n/e is the only physical scale 
surviving the continuum limit. Expanding the rhs of the first line of equation ( |3.13|) at 
small e up to order 4, we get 

1 ( 1 (5 (8»'(0)y 3 flW(O) 

I 32 V^(°) / 32 01(0) 1 1 



to be identified with gr = 1/8(1 — e 4 ) as in eq. (|4.1|) . This fixes implicitly the value of r as 
a function of u by 

4 _ 5 fe'{'(0)\ 2 3 6»f } (0) 



U 32 V^i(0) J 32 01(0) (4 ' 6) 
We may now expand R and R n as given by eq.( |3.13| ) up to order 2 in e to obtain the 
relevant scaling function 

e 2 6"'(0Y 
R = 2 1 1 v ' 



4u; 2 9[(0) J (4.7) 
R n = 2(l- e 2 U(r)) 
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with r as in eq. (|4.1| ) and where the scaling function U{r) is now given by 



U{r) 



4w 2 6>' 1 (0) 



6>i"(0) 




(4.8) 



where we have identified the Weierstrass p function || with half-periods u and uj' = tuj 



where r is fixed by eq.(O), which amounts to 




4 



(4.9) 



3 



where g2 is the first elliptic invariant of the Weierstrass function. The scaling function 
U{r) may be viewed as the potential felt by the random tree in the presence of the two 
walls. 

When sending the second wall to infinity, i.e. taking uj — > oo, we immediately recover 
the above result ( [4.4j ) as g2 = 4/3 fixes the second half-period uj' = z7r/2, in which case the 
scaling function degenerates into ([4.4j ). Taking L = L(g), i.e. uj = 7r/2, fixes uj' = ioo which 
corresponds to q = 0, in which case the scaling function degenerates into a trigonometric 
function 

U( r ) = -4r-,- 1 ( 41 °) 

sin (r) 

The values of uj E [7r/2, oo] (L > L(g)) correspond to taking q £ [0,1], while those in 
[0, 7r/2] (L < L(g)) are obtained for q G 0], where q* = — lim^^oo q c (L) = e _7rv ^. 

Note that we could have derived the scaling limit of R n without using the explicit 
solution ( p.!3| ) by plugging the ansatz R n = 2(1 — e 2 U(r)) in the original equation ( |2~T| ) 
and expanding up to order 4 in e to obtain a differential equation for U, namely 



and require that U{r) diverge at r = and r = 2uj, with no divergence in-between. This 
equation is to be compared with that satisfied by p, namely p" = 6p 2 — #2/2, which allows 
to identify U(r) = Sp(r) provided g% = 4/3. 

5. A simple application: escape probability from a fixed domain for a spreading 
population 

As mentioned in the introduction, rooted planar trees may be used to model discrete 
branching evolution processes. Let us assume for instance that an initial parent individual 



U" = 2(U 2 - 1) 



(4.11) 
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(materialized by the root vertex) gives rise in one generation to a number k > of children 
individuals with probability pk, each child itself independently giving rise to subsequent 
generations with the same probabilities. The resulting genealogical tree is nothing but a 
planar tree (without labels). In the following, we concentrate on the most natural choice 
Pk = (1 — p)p k , with < p < 1, and where the prefactor (1 — p) ensures the correct 
normalization and may be interpreted as the probability of death without descendents. 
This choice is expected to capture all possible physics of the problem, as the parameter p 
allows to explore all possible values of the average number of children p/(l — p), known to 
be the only relevant quantity in the problem |]]. 

We may now turn the branching process into a spatial branching process by allowing 
the individuals to spread in a one-dimensional target space. More precisely, we consider a 
discrete version of the problem in which the individuals may occupy integer-valued posi- 
tions, with the diffusion rule that each child lives at a position differing by ±1 from that of 
its parent, with an equal probability 1/2. Using these positions as vertex labels, we may 
write the following master equation for the extinction probability of a family 



where E n (T) stands for the probability that an individual sitting at position n has no more 
descendent at generation T, with the initial condition E n (0) = 0. Eq.(|0]) is obtained by 
enumerating all possible configurations of the first generation children, and noting that 
the joint probability of extinction of all their descendents is the product of individual 
extinction probabilities before generation T — 1. Comparing eq. (|5.1|) to eq . (pTT|) , we see 
that E n = limy^oo E n (T) obeys the same equation ( |2.1| ) as (1 — p)R n , provided we set 



This allows to identify E n = (1 —p)R n by noting that this choice corresponds precisely to 
the stable fixed point of the recursion relation ( |5.1| ). We may therefore interpret (1 —p)R n 
as the probability of extinction of a family whose first generation's parent sits at position 
n. Accordingly, we interpret 1 — (1 — p)R n as the survival probability for such a family. 

For unconstrained positions (no wall in the former language), we simply have a trans- 
lation invariant probability of survival 



E n (T) 



1 — p 




1- §(£ n+1 (T-l) + £ n _x (T-l)) 




(5.3) 
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obtained by substituting g = p(l —p)/2 into eq. fl2.2p . This result is totally insensitive to 
the diffusion process, and is the same as that for unlabeled trees. It displays a first order 
singularity (discontinuous derivative) at p = p c = -|. This classical result in the theory 
of branching processes |J is a particular case of a more general statement for so-called 
Galton- Watson processes that the genealogical tree is almost surely finite (here S(p) = 0) 
if and only if the average number of children is less or equal to one (here corresponding to 
J2jPj = — P) — 1' namely p < 1/2). 

Apart from the extinction probability E n (T), one natural quantity to study is the 
probability C n (T) for the population spreading from a position n to remain confined within 
a given connected domain T> until generation T. This quantity is readily seen to obey the 
same equation ( |5.1|) as E n (T) for all n in T>, but with a different initial value at T = 0, 
namely C n (0) = 1, and also the condition that C n (T) = for all n outside of V. We 
distinguish the two possible cases where (i) V is a half-line, say [0, oo) or (ii) V is a 
segment, say [0, L]. The only relevant conditions are at the boundary of the domain and 
read respectively (i) C_i(T) = and (ii) C_i(T) = Cl+i(T) = 0. In the limit when 
T — * oo, we may therefore identify C n (oo) = (1 — p)R n with R n given by our (i) one- wall 
and (ii) two-wall solutions of Sects. 2 and 3. The results will be best expressed in terms of 
the quantity S n = 1 — C n (T = oo) which is nothing but the probability for the population 
to escape from the domain. 

In the one-wall case, the solution (|2.3|) leads to the population's escape probability 



l-|2p-l|(l-x- +1 )(l-x"+ 5 ) 

2p (l-x-+ 2 )(l-x-+ 4 ) K } 

where 

_ l-|2p-l|* 
x = x{p) = —j- (5.5) 

1 + 1 2p — 1 1 2 

Note that for any fixed n > the escape probability S n is strictly positive as soon as p > 0, 
and moreover that it still displays a singularity at p = p c = 1/2 but of weaker third order 
type (discontinuity of the third derivative) as is readily seen by expanding S n up to order 
3 in powers of \2p — 1|. Note also that lhm^oo S n = S(p) as in eq.( |5.3| ), expressing the 
equivalence in probability between surviving forever and reaching infinitely distant points. 
Note finally the following simple expression for the escape probability S n at the transition 
point p = \: 

S '^=\)= („ + 2) 3 (n + 4) ' n> -° (5 ' 6) 
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The scaling limit (|4.1| ) may be used to study the vicinity of the transition point by setting 
p = p c (l + rje 2 ) with 77 = ±1 according to whether we approach the transition from above 
or below. Eq.( f4.3| ) allows to interpret the scaling function U{r) as describing the scaling 
behavior of the escape probability S n ~ S n around p = ~, with 

S n = ^U(ne) + „) = |2p - 1| ( -^-±_ + ^ + {2p _ 1} (5 . 7) 

valid in the scaling region of large n and ne = 0(1). This scaling function displays clearly 
the above-mentioned third order transition with 

S n = 4 + (2p - 1) + -(2p - l) 2 - -^|2p - 1| 3 + 0((2p - l) 4 n 6 ) (5.8) 

77/ O Do 




(a) (b) 

Fig. 5: Escape probabilities S n of eq. (|5.4|) as functions of p in the presence of 



one wall, (a) for p e [0, 1] and from top to bottom n = 0, 1, 2, 3, 4 as well as 
n = 00 in which case we recover the no- wall solution S(p) of eq.( |5.3| ); (b) for 
p in the critical region p ~ | and from top to bottom n = 5, 7, 12, together 
with the expected scaling limits S n as given by eq. ( |5.7[ ) with a proper shift of 
n — > n + 3 ensuring the perfect matching of the curves. 



We have represented in Fig. 5 the escape probabilities S n (p) and their limit S(p) as 
functions of p G [0,1]. We have also blown out the critical region around p = 1/2 to 
compare the exact solution with its scaling limit. 

In the case of two walls, the solution (|3.13| ) leads directly to the probability S n = 
1 — (1 — p)R n of escaping from the interval [0, L\. 
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Fig. 6: Plot of the increase in the escape probability Sn ' (p) — S n (p) from the 
one- wall situation to that with two walls for L = 5 and n = 0, 1,2 (bottom 
to top). 

For illustration, we have displayed in Fig. 6 the increase in the escape probability 
S n L ^ (p) — S n (p) from the one- wall situation to that with two walls for L = 5 and n = 0, 1, 2, 
as functions of p E [0, 1]. The increase is maximal at p = ^. 

For finite L, the critical value g c (L) > 1/8 is never attained as g = p(l—p)/2 e [0, 1/8] 
for p G [0, 1], hence the singularity at p = \ is suppressed. However it is restored in the 
scaling limit where L — > oo with L\2p — lj 1 / 2 fixed, as g c (L) — > 1/8 in this case. The 
scaling function U{r) given by eq.(^.8|) again describes the scaling behavior of the escape 
probability S n ~ S n in the vicinity of p = |, with the result: 

S n = 3|2p - l\p{n\2p - I] 1 / 2 ) + [2p - 1) (5.9) 

where the Weierstrass p function must be taken with fixed half-periods 10 = L\2p — l| 1 / 2 /2 
and uj' = tu, such that g2(u,u') = |. We again note that S n displays a third order 
singularity at p = \ by expanding 

S n = A + ( 2p_l) + A^ 2 ( U ;, w ')n 2 (2p-l) 2 + 7 |^3(^,^)n 4 |2p-l| 3 +0((2p-l) 4 n 6 ) (5.10) 

where g<i and gs stand for the elliptic invariants of the Weierstrass function and with 
the constraint that g 2 (ui,u!') = 4/3 which fixes u' and consequently gs as functions of 
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uj = L\2p — l| 1 / 2 /2. Note that the singularity of <S n at p = 1/2 disappears at the modular 
invariant point r = i i.e. q = e _27r , where gs(u>, iu) = 0, which causes all odd powers of 
1 2p — 1 1 to vanish in the series expansion ( |5.10| ) . 

Of particular simplicity is the case when q = in eq. (|3.13| ), corresponding to g = g^ 



as in eq. (|3~TTD , i.e. p = (1 - tan 2 (7r/(L + 6)))/2 or p = (1 + tan 2 (7r/(L + 6)))/2. The 



formula for the associated escape probability S n reads 

s = J ^r^y (»^r-»l)^(-^) - 2sin2 (^)) p=K 1_tan2 (^)) , 511 > 
{ sin (^) sin(7r£ ±|) p- 2 ^+ tan u +6 ;; 

This allows for framing the exact value of S n at p = 1/2 between these two values for 
all L. For large L, we may identify w = lim^^oo L\2p — l| 1 / 2 /2 = n/2 for both values 
of p = (1 ± tan 2 (7r/(L + 6)))/2, and r = zoo to ensure g = 0, in which case the scaling 
function reduces to Q4.10|) and therefore eq. flOl) turns into 

S " = |2P - 1| U W 2 P 3 - 1 |^) - 1 ) +(2P - 1) (5 ' 12) 



in agreement with the large L limit of eq.(|5.11|). 



6. Another solvable case of tree embedding: dilute SOS model on a random 
tree 

We may consider a slightly different version of labeled trees, in which we impose the 
weaker constraint that any two adjacent vertices of the tree must have labels differing 
by ±1 or 0. As shown in Ref. ||, this version of labeled trees is that involved in the 
enumeration of tetravalent planar graphs. In the language of spreading of a population, a 
child may now stay at the same position as its parent. This corresponds to a dilute SOS 
version of the case studied in this paper. 

The main recursion relation is now replaced by 

R n = ; r (6.1) 

1 — g{Rn+l + Rn + Rn-l) 

where R n is the generating function for rooted trees with root vertex labeled by n e Z, 
and a weight g per edge. Again, we may consider three types of boundaries: no wall, one 
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wall and two walls. The no- wall case is easily solved, with all the R n 's equal to the solution 
of R = 1/(1 - 3gR) with R = 1 + O(g), namely 



R = R ig) = 1 -^IEm ( 6 . 2) 

The one- wall case corresponds to setting = and only considering the i? n 's for n > 0. 
It was solved in Ref. H , with the result 

R n = R , u n = x 2 —x 2 (6.3) 

with = R(g) of eq.(|6.2|), and where x is the solution oix + l/x + A = l/(gR) with, say, 
modulus less than 1. Note the slight difference in the index shifts when compared with eq. 
). The main recursion relation reduces this time to a quartic equation for the it n 's: 

u n u n+1 u n+2 u n+3 = ^u 2 n+1 u 2 n+2 + gR(u n _iul +2 u n+3 + u 2 n u 2 n+2> + u n u 2 n+1 u n+A ) (6.4) 



supplemented by the initial condition u-i = 0. Eq.(|6.4j) is easily checked for all k by 
setting 1/R= (x 2 + x + l)/(x 2 +4x + 1) and gR = x/(x 2 + 4x + 1). Finally, in the two-wall 
case where we require R-i = Rl+x = and only consider n = 0, 1, 2, L, we have found 
the elliptic solution 

u n u n+3 

tin — ti~ 



Un+lUn+2 (6.5) 

u n = 0i((n + l)a) 

with x = e 2l7ra and 6*i as in eq. (p^) , and which solves eq. (|6.4| ) provided we take 

R _ A e 1 (a)e 1 (2a) f 61(a) 1 9[(2a) 



'^(0)0! (3a) \^i(a) 2^i (2a). 

6'i(0) 2 6' 1 (3a) ( 6 -6) 



16^i(a) 2 ^i(2a) 



L (a) 2 ^ (2a) 



The boundary conditions are again satisfied for two choices of the parameter x: (i) x = 
e 2i7r/(L+5) anc | ^jj-j _ ^i/(l+5) ( w h en q > q), the latter leading to the same physical 

solution as the former by modular invariance. Picking again the first solution, we must 
take a = 1/(L + 5), and may view the equations for the solution as parametrized by q. 
This leads to the continuum limit, upon taking g = (1 — e 4 )/12, n = r/e, 2u> = (L + 5)e 
and R n = 2(1 — e 2 U{r)). We end up with a scaling function U(r) = 2p(r) in terms of 
the Weierstrass p function with half-periods uj as above and uj' fixed by now requiring 
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that g2(uj, uj') = 3. We also recover the one- wall case by taking the limit L — > oo, namely 
uj = oo while u/ = z7r/\/6, in which case U{r) = 1 +3/sinh 2 (^372r), a result already 
obtained in Ref. ||. 

Note finally that all scaling functions coincide with those of Sect. 4 up to a global rescal- 
ing r — ► yfr and — > v/f^ - This confirms the expected universality of the continuum 
limit. 



7. Conclusion 

In this paper, we have extensively studied a model of random rooted planar trees em- 
bedded in a discrete one-dimensional target space. In particular, we have derived explicit 
expressions for the partition function of the model with various target spaces, namely the 
whole integer line, a half-line, and a segment. To obtain these, we have shown that the 
partition functions actually obey recursion relations and that the particular target at hand 
translates into various boundary conditions. We have also derived the corresponding scal- 
ing functions in the continuum limit for which the recursion relations turn into differential 
equations. A different approach, popular among probabilists, consists in studying directly 
the continuum limit of embedded random trees in the form of continuum spatial branch- 
ing processes 0, giving rise to partial differential equations. It should be possible in this 
context to solve these equations with wall-type boundary conditions, and to recover our 
continuum results. Some work in this direction appeared recently |10|, where an analogue 
of our one-wall case can be found. 

The striking simplicity of the solutions (|2.3|) (|3J]) as well as (|6.3|) dT5|) are directly 
linked to the "integrability" of the corresponding non-linear recursion relations ( [2.1| ) and 
respectively. One possible explanation for the integrability uses the interpretation of 
labeled trees in the context of planar graph enumeration. More precisely, as shown in detail 
in appendix A below, the SOS and dilute SOS models on trees respectively occur in the 
enumeration of rooted planar Eulerian triangulations (i.e. triangulations with bicolored 
faces) and of rooted planar quadr angulations. In both cases, the labels n correspond to 
geodesic distances along the graph from the root. This reformulation suggests a possible 
connection with matrix models, known to be integrable. As already hinted in Ref. @ in 
the context of graph enumeration, the equations ( |2.1[ ) and ( |6.1| ) are very similar to those 
arising in the context of the matrix models used for generating Eulerian triangulations 
and quadrangulations respectively, namely R n = (n/N)/(l — g(R n+ i + R n -i)) and R n = 
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(n/N)/(l — g(R n+ i + R n + i? n _i)), where N is the size of the matrices. The remarkable 
point is that the index n is no longer related to the geodesic distance along the graphs, but 
to their genus. Soliton theory seems to indicate that integrability survives when changing 
the recursion into R n = (a+(3n)/(l— g(R n+ i+R n -i)) or R n = (a+(3n)/(l—g(R n+ i+R n + 
R n -i)) for some constants a and /?, which could interpolate between the two problems. 

Beyond the two (quadratic) examples of this paper, we have at our disposal a host of 
non-linear recursion relations all used for enumerating possibly decorated planar graphs 
while keeping track of some geodesic distances, and which were found to be integrable 
as well (see Ref.[[| for details). The solutions display some multicritical behavior corre- 
sponding to higher order critical points, with non-trivial hierarchies of scaling functions. 
It would be interesting to find some proper interpretation of those equations in the context 
of embedded trees or alternatively of population-spreading processes. 

It would also be interesting to classify the target spaces leading to integrable models 
of embedded trees, or alternatively to spot among all possible discrete spatial branching 
processes those with integrability properties. We may then hope, by changing the nature 
of the discrete target space, to be able to reach new critical points. 

Appendix A. Graph interpretation 

As mentioned above, well-labeled trees, i.e. trees with non-negative labels correspond- 
ing to the one- wall situation in the dilute SOS version of Sect. 6 were introduced in Ref.@ 
in the context of graph enumeration. More precisely, it was shown that there exists a 
bijection between these well-labeled rooted planar trees and rooted planar quadrangula- 
tions. This allows to interpret the quantity R n of eq.( |6.3| ) as the generating function for 
planar quadrangulations with both a marked (origin) vertex and a marked oriented edge 
linking a vertex with geodesic distance m from the origin to a vertex with geodesic dis- 
tance m + 1 from the origin with m < n, and with an activity g per vertex. Beside this 
equivalence, there exists yet another bijection, now in the dual language, between rooted 
planar tetravalent graphs (dual to the above quadrangulations) and decorated (so-called 
blossom-) binary trees 0. This bijection was extended so as to keep track of geodesic 
distances between faces in Ref. [RJ . In this language, the generating function Z n for two- leg 
planar tetravalent graphs with the two legs distant by at most n was shown to obey the 
recursion relation Z n = 1 + gZ n (Z n+ i + Z n + Z n _i), with Z-\ = 0, a direct consequence 
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of the above bijection with blossom binary trees. This equation is nothing but yet another 
form of eq.( |6.1|) and allows to identify Z n = R n of eq.(|673|). 

Remarkably, our slightly simpler equation ( |2.1| ) also admits two analogous interpre- 
tations in terms of graphs, now related to the enumeration of rooted planar Eulerian 
triangulations (i.e. triangulations with bi-colored faces, say in black and white) or dually 
to that of rooted trivalent bipartite planar graphs (say with black and white vertices). In 
this dual language, we may rely on a bijection [|| between rooted trivalent bipartite planar 
graphs and properly decorated binary trees. Keeping track of the graph-geodesic distances 
in these binary trees leads to the equation Z n = 1 + gZ n (Z n+ i + Z n _i) for the generating 
function Z n of two-leg trivalent bi-colored planar graphs with the two legs attached to 
vertices of opposite colors and at geodesic distance at most n. In this approach, the proper 
definition of the geodesic distance on the graphs makes use of oriented paths linking faces, 
with the constraint that a step across an edge between two faces always leaves the black 
vertex on the left. The equation for Z n is yet another form of eq.( |2.1| ) which allows to 
identify Z n = R n of eq.(|2.3|). 

In terms of triangulations, a bijection similar to that of Ref. || may be established be- 
tween rooted planar Eulerian triangulations and the well-labeled (SOS) trees corresponding 
to the one- wall situation of Sect. 2 as follows. First we replace the face-bicoloration by the 
compatible orientation of all edges in such a way that each triangle is either clockwise- 
or counterclockwise oriented. This allows to define the geodesic distance from a vertex to 
another as the length of a minimal path respecting the orientation of edges. Picking as 
origin some vertex, a well-labeled tree is obtained by retaining for each clockwise-oriented 
triangle the edge linking the two farthest vertices from the origin, and labeling each vertex 
by its distance from the origin minus one, as illustrated in Fig. 7 (see Ref.|J for details 
and proofs). This allows to interpret R n as the generating function for planar Eulerian 
triangulations with a marked (origin) vertex and a marked oriented edge linking a vertex 
with geodesic distance m from the origin to a vertex with geodesic distance m + 1 from 
the origin with m < n, and with an activity g per vertex. 

In the language of graphs, we may use our two- wall solutions to enumerate bounded 
graphs as follows. The quantity 

cf^^-e^ (a.i) 

where Rn is the solution of eq. (|2.1|) (resp. ( |6.1|) ) with two walls at positions —1 and L + l, 
is the generating function for Eulerian triangulations (resp. for quadrangulations) with a 
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(a) (b) 

Fig. 7: A sample Eulerian triangulation (a) with a marked origin vertex (la- 
beled 0) and a marked oriented edge (big empty arrow). We have indicated 
for each vertex its geodesic distance from the origin (respecting the edge- 
orientations). A labeled rooted tree (b) is obtained by retaining for each 
clockwise-oriented triangle the edge connecting the two farthest vertices from 
the origin, and picking as the root the end of the previously marked edge. 
The vertex labels on the tree are simply the distances of the graph vertices 
from the origin minus one. 

marked (origin) vertex, a marked oriented edge linking a vertex with geodesic distance n 
from the origin to a vertex with geodesic distance n + 1 from the origin and which are 
"bounded" in the sense that the geodesic distance of all the vertices from the origin is less 
or equal to L + 1. As an example, the rooted Eulerian triangulation of Fig. 7 contributes 
to G[ L) for all L > 2. 
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